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A generalization of the Zerilli master variable for a dynamical spherical spacetime 
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The evolution of polar perturbations on a spherical background spacetime is analyzed. The 
matter content is assumed to be a massless scalar field. This provides a nontrivial dynamics to 
the background and the linearized equations of motion become much more involved than in the 
vacuum case. The analysis is performed in a Hamiltonian framework, which makes explicit the 
dynamical role of each of the variables. After performing a number of canonical transformations, 
it is possible to completely decouple the different perturbative degrees of freedom into constrained, 
pure-gauge and gauge-invariant variables. In particular, two master variables are obtained: one 
corresponding to the polar mode of the gravitational wave, whereas the other encodes the complete 
physical information about the perturbative matter degree of freedom. The evolution equations for 
these master variables are obtained and simplified. 


I. INTRODUCTION 

The study of linearized perturbations of known background solutions of general relativity has had very important 
contributions to our current understanding of different gravitational scenario. As in other field theories with con¬ 
straints, one of the main problems of this approach is the identification of physical degrees of freedom. There are 
usually two approaches one can follow for such a purpose. On the one hand, it is possible to impose convenient 
gauge-fixing conditions and work on a certain gauge. On the other hand, one can construct gauge-invariant variables 
so that any physical result is unambiguous and valid in any gauge. 

The Hamiltonian formalism of general relativity gives a very clear and transparent notion of the gauge dependence. 
In particular, the Hamiltonian is a linear combination of first-class constraints, which are the generators of gauge 
transformations. In the context of perturbation theory, such a Hamiltonian formalism was pioneered by Moncrief [l| 
to study the nonspherical perturbations of the Schwarzschild black hole. In that reference it was shown that, for this 
solution, it is possible to perform several canonical transformations explicitly in such a way that the initial twelve 
perturbative variables (the six components of the perturbed spatial metric in combination with their corresponding 
conjugate momenta) are reorganized into two physical pairs, which encode the complete physical information of the 
gravitational wave, and four gauge pairs. In each gauge pair, one of the variables is constrained to vanish on shell, 
whereas its conjugate variable is nonphysical. Furthermore, the two physical pairs obey unconstrained evolution 
equations and are equivalent to the Regge-Wheeler Q and Zerilli Q master variables. In this way, the physical 
degrees of freedom are explicitly decoupled from the gauge degrees of freedom and the dynamical behavior of the 
system is completely described by the two physical pairs. This technique was also applied to other specific solutions 
of Einstein equations like Reissner-Nordstrdm Si: Oppenheimer-Snyder S or Friedmann-Robertson-Walker Si- 

Regarding spherically symmetric background metrics, it is well known that the perturbations can be classified into 
two different sectors (axial and polar) depending their polarity. At linear order these two sectors decouple, but at 
second and higher orders they interact as, for instance, the coupling between two axial modes give rise to both axial 
and polar modes i . For the Schwarzschild metric the Regge-Wheeler S and Zerilli Q variables encode respectively 
the axial and polar physical degrees of freedom. Both are master scalars since they obey unconstrained evolution 
equations and the complete perturbed metric can be reconstructed in terms of them. These two variables were initially 
defined on a fixed perturbative gauge but, as commented above, were later obtained by Moncrief on a generic gauge. 

A very convenient framework to study perturbations around generic (possibly dynamical) spherically symmetric 
spacetimes was presented by Gerlach and Sengupta [iS[ni- This is a very geometrical framework, where the four 
dimensional manifold is decomposed as the product between a two dimensional Lorentzian manifold with boundary 
and the unit two-sphere. The construction of perturbative gauge-invariant variables is explicitly performed and, for 
the axial case, a master scalar variable is constructed. This master scalar obeys an unconstrained wave equation and 
can be coupled to any kind of matter. Therefore, it can be considered as the generalization of the Regge-Wheeler 
variable to dynamical spacetimes. On the contrary, for the polar sector, there is no known master variable valid for 
any spherically symmetric background. Nonetheless, some particular results for specific metrics have been obtained. 
For instance, on a vacuum background, the gauge-invariant combinations of the linearized stress-energy tensor have 
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also been included in On the other hand, in Ref. 0 a polar master scalar was dehned, for a vacuum 

background, which was later generalized to nonlinear electrodynamics [Tsjl for any background solution. 

The present paper is motivated by the search of a unique gauge-invariant master variable encoding the polar 
gravitational wave for any matter model. In principle it is not clear that it can be constructed but, apart from its 
conceptual relevance, such a polar master variable would be of great use for several applications. For instance the 
numerical resolution of an unconstrained equation is, in principle, much easier and precise than the resolution of a 
coupled system of several equations which, in addition, are subjected to certain constraints. In addition, it could 
be particularly useful in the very hard problem of matching of polar gravitational waves through moving surfaces, 
for instance the surface of a supernova explosion [l6l - [l8l| . This problem becomes even harder beyond first-order 
perturbation theory Il9l-l22ll . 

In a previous paper [23[ the Gerlach-Sengupta axial master scalar was reobtained , making use of the Hamiltonian 
techniques explained above, for a dynamical spherical background spacetime with a matter content of a scalar held. 
Here the analysis of that paper is reproduced for the polar sector, as an example of a dynamical spacetime for which 
no polar master variable is known. 

The rest of the paper is organized as follows. In Sec. |H]the Hamiltonian framework for linearized perturbations on a 
generic background is briehy reviewed. Section Hill introduces the notation and the equations of motion corresponding 
to the particular background we will be dealing with; a spherically symmetric spacetime with a massless scalar held. 
In Sec. llVI the polar part of the perturbative variables are decomposed into tensor spherical harmonic and a number 
of canonical transformations are performed in order to decouple the gauge and the physical degrees of freedom. In 
Sec. |V] the evolution equation for the master variables are presented and simplihed. Finally, Sec. |Vl] discusses the 
main conclusions. 


II. GENERAL RELATIVISTIC PERTURBATION THEORY ON A HAMILTONIAN FRAMEWORK 


A. Hamiltonian framework for general relativity 


Let us assume general relativity coupled with a massless scalar held $. In order to perform a Hamiltonian analysis 
of this system, the usual 3-1-1 decomposition of the spacetime is performed. Greek indices will be used for four¬ 
dimensional objects and Latin indices for three-dimensional ones. It is possible to choose coordinates adapted 

to the foliation so that the three-dimensional metric dehned by projecting the four-dimensional one to the 

spatial slices, is given as 


9ij ^ ^9ij- 

Furthermore, the lapse function a and the shift vector j3i are dehned as follows. 


-2 (4) tt 

a := 'g , 


3 . := 


The action for the system under consideration is then given by, 

S= [ dt [ (fx + n$.t -an- , 


( 1 ) 

( 2 ) 

( 3 ) 


where H*-^, which is related to the extrinsic curvature, is the conjugate momentum of the spatial metric gij whereas 
H is the conjugate momentum of the scalar held. As it is well known, the lapse and the shift are Lagrange multipliers 
associated to the Hamiltonian n and momentum constraint ni respectively. These take the following form in terms 
of the basic variables: 


1 


n = 

dg 

n^ = -2Djm 


- - {n^iY 




1 /H 


fJ.g 




n$„ 


where fig := sj det gij and Dj is the covariant derivative associated to gij. 


( 4 ) 

( 5 ) 


B. Linearized perturbations 

Let us begin by dehning a one-parameter family of spacetimes {A4{£),g^,y{e)), where £ is a dimensionless parameter. 
The e = 0 member of the family is referred as the background spacetime. The idea behind linear perturbation theory 
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is to perform a linearization around a background, which is a known exact solution of the Einstein equations. This 
is done by performing a Taylor expansion on the parameter e of all quantities that appear in the equations. Then, 
terms higher than second order in e are dropped. For convenience, we define the operator 


d^F{e) 


(FF{e) 

de” 


£=0 


and, making use of it, introduce the following notation for the perturbative variables: 


( 6 ) 


C:=5a, 

As it was shown in Ref. ESI , the second variation of the action. 



dx‘^ 


+ sn - C6{n) - 




(7) 

( 8 ) 


provides an action functional for the linearized perturbations. (The explicit form of the first and second variations 
of the Hamiltonian and momentum constraints can be found in (^.1 That is, the variation of this last action with 
respect to different perturbative variables 0 gives the linearized Einstein equations. In particular, the variation 
with respect to the perturbation of the lapse C and the shift leads to the constraints obeyed by these linearized 
variables, 


d(H) = 0, 5{n^) = 0. (9) 

These are first-class constraints and, therefore, the generators of the perturbative gauge transformations. The idea of 
this paper is to perform canonical transformations of the perturbative variables so that each of these four constraints 
is simply expressed as one of the new variables. In this way, the different dynamical sectors would get decoupled. 
From the eighteen variables under consideration [the sixteen that appear in ([7]) , in combination with the two functions 
((5$,dn) encoding the scalar degree of freedom], four of them would be constrained to vanish. The four canonical 
conjugate variables of the constrained ones will be pure gauge, whereas the four functions (C, 5*) are Lagrange 
multipliers with vanishing conjugate momentum and thus non-dynamical. Finally, the remaining six variables will 
be automatically gauge invariant. These latter six variables stand for the three physical degrees of freedom of this 
problem: two corresponding to the gravitational wave and one to the perturbations of the matter scalar field. This 
separation between the gauge and the physical sectors is not obvious and can only be performed for highly symmetric 
backgrounds. 


III. SPHERICAL BACKGROUND 


We will consider a spherically symmetric background, which can be decomposed as Ad = x 5^, where is 
a two dimensional background with boundary and the unit two-sphere. Arbitrary coordinates (t, p) on and 
spherical coordinates x°‘ = (0, (p) on are chosen. The lower-case Latin indices stand for coordinates on the two- 
sphere. Due to the symmetry, the lapse can only be a function of the coordinates on M^, that is, a = a(t, p); whereas 
the shift vector has vanishing angular components /3* = (/3(t, p), 0, 0). In this way, the four-dimensional background 
metric takes the following form: 

(ds^)4 = —a^dt^ + a?{dp + jddtY + r'^dVl^. (10) 


Furthermore the following three variables are defined, which completely encode the information contained in the 
background moments H*-^ and H: 


a'^WP 

Hi := -, n2 := - 



Pg pg Pg 

For completeness, here the symmetry-reduced background constraints are provided. 


^ = n. (A-n,)-».« + l(n/ + U)=o, 

n2 + n3$' = o, 

a fig r 


(11) 


( 12 ) 


(13) 
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where prime stands for the derivative with respect to p divided by the function a: 


/' = 



These constraints generate the following evolution equations: 


- [a,t - (^a),p] 
a 

-{r,t - jir^p) 

a 

a 

— (Ill t — /Sill p) 
a 

-(n2.t-/3n2.p) 

a 

-(n3y-/3n3,p) 

a 


2 (Hi - n 2 ). 

--Hi, 

2 ’ 

n3 


sn? 




1 r' {a^r)' 1 


4 r 

2 , tt 2 


Hi - ihn2) 


K 


n? + $" 


2a'r' 2(ar)'' 


ar 


ar 


n3(ni + n2) {ar'^^'y 

^ 4” o * 


1 

2 




(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


The case studied previously by Moncrief is vacuum in Schwarzschild coordinates, which is recovered by choosing 
$ = Hi = 112 = 113 = 0. This restriction greatly simplifies the above equations of motion and in particular both the 
Hamiltonian (|12l) and momentum constraints m are trivially fulfilled. 


IV. POLAR PERTURBATIONS 


A. Expansion in harmonics 


In order to take advantage of the background spherical symmetry we will make use of the tensor spherical harmonics 
in order to decompose the perturbations. Properties and precise definitions of the harmonics that will be used here 
can be found in Ref. . Following their behavior under a parity transformation, different tensor spherical harmonics 
can be divided into two groups: axial harmonics, with a polarity (—1)*+^, and polar harmonics, with a polarity (—1)^ 
Since at a linearized level these two groups of harmonics decouple, it is possible to consider the axial and the polar 
problem independently. As commented in the introduction, in Ref. ( 2 ^ the axial case was developed, whereas here 
we will focus on the polar case. The decomposition of the polar part of the different perturbations is given as follows, 


hijdx^dx^ 

— Pijdx^dx^ 
d‘g 

C 


Bidx'' 





00 I 


E 

1^0 

E 

m— — l 

aHH2)T 

Yi^dp^ + 2{hi)TdpZY^adx‘^ 

+ [Kl^^abYr + GY^Zy^ab] dx^dx\ 

( 21 ) 

CO 

E 

1^0 

1 

E 

m— — l 

a^iPH)T 

Yrdp^ + 2{Ph)TdpZradx‘^ 

' + [{PK)TlabYr + {PG)TZrab] dx<^dx\ 

( 22 ) 

CO 

E 

1^0 

1 

E 

m— — l 


m\rm 

1 5 



(23) 

00 

E 

1 

E 

m— — l 


'^dp+ihoTzr 

adx°-, 


(24) 

00 

E 

1^0 

1 

E 

m— — l 

pTYr. 




(25) 

CO 

E 

1=0 

1 

E 

m— — l 

pTYr, 




(26) 


where the notations by Regge-Wheeler and Moncrief have been followed for the different harmonic coefficients. 
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Since all perturbations are decoupled at linear order, from here on, the {l-,rn) labels from the harmonic coefficients 
and the harmonic tensors will be removed. In addition, we define the following shortening for the sum that appears 
in all decompositions above, 

OO I 

l^m l—O m— — l 


B. Effective action 


In order to obtain the effective action in terms of the harmonic coefficients, decomposition (1211126^ is introduced in 
expression ([8]). The angular part can be integrated by making use of the properties of the tensor spherical harmonics 
as shown in appendix of Ref. . In this way, it is easy to obtain the following form for the effective action for the 
polar part of the linearized perturbations: 


1 

2 


(^2^) polar 



- 


polar 


(28) 


dt{Fo[-aHo/2] + F,[H,] + F 2 [ho]} 


where the dots stand for terms coming from the second perturbation of the background constraints, which do not 
enter the gauge transformations, and the functionals {Fq, Fi, F 2 ) will be defined below. The conjugate momenta are 
related to the harmonic coefficients given by the expansions (I21II26I) in the following way, 


Pi = 

2l{l + l) 

I 

(29) 


a 


Pi = 

ar Pfj, 

(30) 

P3 = 

2ar^P;^, 

(31) 

Pi = 

Xar^Pa, 

(32) 

P = 

ar^p*, 

(33) 


where the star stands for complex conjugate and 

1 (^ + 2 )! 

■ 2 (/- 2 )!’ 


has been defined. 

The polar part of the harmonic decomposition of the linearized constraints is given by, 


m] = 

l,m 
P3 


ni(ni-n 2 )- 


P + l + 2 _ 2i_ 

2fXg _ 


- 2H2'- + ^(Ri - n 2 ) + [-uj -ul + 

r ar^ 2 


rRi - 


ar^ 




{i-m + 2) 


^ i^r\i I I P 


2l{l + l) 


IsiHg] 


= Y.P9y\- 


K +^ir‘^K'y + + £yn3 + [ra-^hi)' - 

2 (a"V2)' , Pi , 2p3 r' p . „ 2 


ar^ 

^2t 




^+%-+^ - Ml - 4 (^"ni )'772+^^^^n2(r^ X)'+n3(^'}>, 


ar^ 


5[nar^^^=Y,PgZa 


-1 {Gpi)' P3 


Pi , ~ 1)(^ + 2) 


l{l + l) 


l{l + l) 


n2G + niH2 - 


r 

2 ( r^IIi/ii 


a 


+ n3(p \ . 


With these relations at hand, the three generators of polar gauge transformations can be written as 


Foif] = 

1 dx^fYS[n], 

(34) 

Fi[f] = 

1 dx^fY^d[Hp], 

(35) 

Fi[f] = 


(36) 
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which act on any smooth arbitrary scalar field /. It is possible to calculate the Poisson brackets between different 
generators, 


{Fo[fi],Fo[f2]} = ! dpar^hf^- 

J Cl fig 

{Fo[fi],Fi[h]} = J dpa/i(^rV2^) , 

{i"o[/i],i"2[/2]} = -Kl + 1) f dpahh — , 

J Mg 

m[/i],i"i[/2]} = j dpar\hf'^- 

J CL jig 

{Fi[h],F2[f2]} = 0, 

{F2[fl],F2[f2]} = l{l + l) [ dpa(/i/2-/(/2)-^, 

J a Pg 


(37) 

(38) 

(39) 

(40) 

(41) 

(42) 


all vanishing on-shell, which confirms that they are first-class constraints. 


C. Canonical transformations: gauge-invariant variables 


Moncrief isolated the Zerilli variable after two canonical transformations on the four pairs 
(/ii,pi), {H 2 ,P 2 ), (KjPs), {G,P 4 ). Here the additional pair {<p,p) for the scalar field is also present; and the 
fact that the background is dynamical makes the problem harder. We will instead proceed in five steps, to clarify 
the role of each step and simplify the computations. In particular we will first eliminate the two gauge degrees of 
freedom related to the momentum constraint (which are rather trivial and very similar to the axial case) and then 
remove the gauge degree associated with the Hamiltonian constraint, the nontrivial step of this computation. 

There are many possible transformations that implement this program, but we would like them to obey certain 
minimal criteria. First, they should be algebraic transformations so that they do not involve any integration in the 
process and can be performed explicitly. Second, they should not require dividing by any background object that could 
vanish, in particular one of the background momenta (ni,n 2 ,n 3 ). And third, in order to obtain a generalization of the 
Zerilli master variable, all transformations should be well defined in the vacuum limit, that is, when taking vanishing 
values for the scalar field and its perturbations. The full transformation that will be proposed here completely fulfills 
the first and third criterion, but it is unclear whether it also satisfies the second one, as will explained below. 

The first canonical transformation is motivated by the Gerlach and Sengupta choice of gauge invariants 


ki 

k2 

ks 

ki 


/C5 


K + 


l{l + l) 
2 




hi 


G, 

„2 

a-^hi - -G', 

if - 



(43) 

(44) 

(45) 

(46) 

(47) 


TTl 

= P3, 

(48) 

7r2 

= P2, 

(49) 

TTS 

l{l + 1) 1 /• 2 \l 

= P4 ^ P3 ^a{r pi) , 

(50) 

714 

2r' 

= api - 2a(a“ V 2 )' H- P3 + P 4>', 

r 

(51) 

T^5 

= p- 

(52) 


which requires the canonical momenta 
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In terms of these new variables, the components of the perturbed momentum constraints are written as 


1 








— Illfcn ~ 2^2 




a fXg 


+ 




n 2 {k4 - ksrr') + fc4n3$" + k 4 ^{ry + HgA:'| 


] polar 




27ro 

+ felli + Ilafcs -5~n2fc4 


1(1 + 1) 


- \{Tl4r%y + 

/(( + !) a 


(53) 
• (54) 


The explicit form of the perturbation of the Hamiltonian constraint <5 [77] in terms of these variables is not displayed 
because it is a very lengthy expression and does not contribute in any way to the present discussion. 

A second canonical transformation is performed, which converts the momentum constraints into canonical variables. 
Because of the requirements we want to impose in all our transformations, only the momenta 744 and 743 can replace 
the constraints (15^ and (IMl) respectively. 


741 — ar^ (^2fc4)^ 

(55) 

772 + + ar^Hifc^ — a(r^ni)'fc4. 

(56) 


(57) 


(58) 

775 + Eillai.2n3fc3 - a(r^n3A:4)'. 

(59) 


74l = 

7^2 = 

^3 = 

744 = 

^5 = 

The division by the tensor harmonics must be understood just as removing them, as well as the summation symbol, 
from the above expressions (l53lj54ll . These last transformations for the momenta do not affect the position variables, 

(60) 

(61) 

(62) 

(63) 

(64) 


ki = ki, 

k2 = k2, 
h = ks, 

/t4 5 

k^ — ^5 ■ 

In terms of these last variables, the perturbative constraints take the following simpler form, 

2 ^-2^/ 


m] = -n 


’"I + (Hi - n2)^ + n3^ + - ^^2 + ~kk^' 


+ 


ni(n2-ni)- 


(l-l)(l + 2) 


-n 


fci - 


ni(n2 — Hi) + 


(1 — i)(i + 2) + 4 


42 f , 


-s[ng] 

a 


ma 


I polar 


l,m 

l,m 


773 


l{l + l)ar‘^' 


(65) 

( 66 ) 

(67) 


The gauge freedom contained in the perturbed momentum constraints has been fully isolated. The variables k^ and 
k 4 are gauge dependent and non-dynamical because their conjugate momenta 743 and 744 are constrained to vanish. 
We are left with a system of three degrees of freedom (fci, k 2 , fcs) and the single constraint (IS51) . 

Following the same procedure, at this point one should make another canonical transformation and convert the 
Hamiltonian constraint into one of the variables. Because of the first criterion we want to impose, we can not convert 
any of the variables {fci, k 2 , ^ 5 } which appear in (I65p into the full constraint. But, because of the second requirement, 
we can neither do it for any of the momenta { 741 , 742 , 745 }. Therefore, the idea is to first make a transformation that 




















removes second-order derivatives of ki from the constraint dMl), so that all first derivatives of the perturbed objects 
can be absorbed in a single term. We use an arbitrary constant 7 to parameterize the transformation, 


which will introduce a first derivative of 1:2 


TTl 

^2 

^3 

7f4 

^5 


ki = fci, 

h h h 

k 2 = ^2-;-7^1? 

r 


( 68 ) 

(69) 

(70) 


^3 = ^3, 

^4 = ^ 4 , (71) 

h = h, (72) 

in the Hamiltonian constraint through the transformations of the momenta, 


= TTl -I- (7 - l) 7 r 2 - I 


T^2 


= 7r2, 
= ^3, 
= 7r4, 
= ^5- 


(73) 

(74) 

(75) 

(76) 

(77) 


In this way, the Hamiltonian constraint can be written as a sum of a full derivative term (which will be later 
promoted to the polar gauge-invariant geometric master variable) and a linear combination of variables ki and TTi with 
no derivatives. 


m] = 


l,m 


•7f9 ~ 1 ~ 1 ~ 

- 2^fc2 -b - -iV^ki 

av^r r r' 


-Hi 


TTl 

9 

ar^ 


H. 


7^3 

> 2 
ar^ 


_ Hi - n2 + yHi) ^ 

r / ar^ 


T T 

-7 + ^ 7 ^^ 


x 

- - ( 1 / 


3(r')2 + ^)ni(n 2 - Hi) - 


l{l + 1 ) , , 2 

^2 (^ + 0 ') + ^(^~ 0 ') 


/ci 


3(r')^ 


42 ^ , 


(78) 


where the background potentials, 

^:= ^7 +ni(n2-ni) + —, := h +(27-3)^^, (79) 

have been defined. As it has been anticipated, this clearly motivates another canonical transformation in which the 
term that appears inside the full derivative [in the first line of Eq. (l78ll ] replaces the canonical variable k 2 - Note that 
using ki instead to replace such term would require dividing by Hy, which is a background object that could vanish, 
whereas using k^ would not provide a well defined vacuum limit. The fourth canonical transformation takes thus the 
following form: 


ki 

= ki, 

(80) 

k2 

7 f 2 r' ~ 1 1 ~ 

= Hi 2 fc2-f $ fcs .Hy/ci, 

ar r r r r 

(81) 

h 

= h, 

(82) 

ki 

= ki, 

(83) 

h 

~ ^5; 

(84) 

^1 

II 

1 

to 

to 

to 

(85) 

772 

rTT2 

( 86 ) 

^3 

= -^3, 

(87) 

71-4 

= "^4, 

( 88 ) 

^5 

= ^5 + ^^2. 

(89) 
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Now the Hamiltonian constraint does not contain 1:2 and, in addition, it has no explicit dependence on the constant 
7 . But, more importantly, we have achieved what we were looking for: it neither contains derivatives of fci. 


(5[iW] = V MgF <{ - Hi^ + Ha^ + 4 + 


/ , I 
Z,; 




•^5 


h 

2 


V 


+ 3r' I - k 


‘J>" + y(7 + 3-' 


(90) 


where we have defined the background coefficient 


B = 


(r^vy 


rr 


2(r'y 




1(1 + 1 ) 


+ ni(n2-ni) + <^)i? 


2^ii±ll + nl 


(91) 


This fact permits us to perform the final fifth canonical transformation, which converts the Hamiltonian constraint 
into the first of the variables of the problem. 


Qi 

Z 

Qs 

Qa 


Pi 


Pz 

P 3 

Pa 


Dh - + Ha^ + 4 + ^ ( ^ + 3r' ) - 

ar^ ar^ 2 

k2, 

h, 
h, 


TTS . T, k2 fV 


fcs + 

D' 

^2 + 

^3, 

■^4) 

^5 


H; 


3 _ 

TZ"!, 


ar^D 




D 


fv 




(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

( 100 ) 

( 101 ) 


At this point we have succeeded in separating the physical degrees of freedom (Z, Pz) and (cj), P^) from the gauge 
degrees of freedom (Qiyi), (Qs^Ps) and (Qi,Py). However in this last transformation the background object D 
appears as denominator and it is not clear to us whether this object can vanish or not. In vacuum, for a Schwarzschild 
solution, defining A := (I — 1)(Z + 2)/2, we have 


D = 


1 1(1 + 1 ) 

1 - 2M/r r3 


(Ar + 3M), 


( 102 ) 


which is always positive. It is reasonable to assume that for spacetimes close enough to Schwarzschild (though possibly 
dynamical), the variable D will also be positive. If this was the case we would have succeeded in implementing to 
completion the procedure while obeying the three imposed criteria. If not, analyzing the procedure that has been 
followed, it seems quite difficult to achieve the construction of gauge-invariant master variables without dividing by 
a never vanishing background object for a generic background gauge. Note that, once the canonical transformation 
(I431I52I1 is performed, there is no much freedom left in the procedure if one insists on imposing the three criteria: 
performing only algebraic transformations, not dividing by a possibly vanishing background object and having a well 
defined vacuum limit. More specifically, the momenta tts and 7 r 4 are the only variables that can be used to solve 
for the constraints (l54l) and (l5^ respectively. This leads to the form (l 6 ^ of the linearised Hamiltonian constraint. 
In that expression none of the variables can be used to solve the constraint algebraically. Thus, next parametrized 
transformation (j68H77l) is performed in order to concentrate all derivatives in a unique full derivative. The term inside 
this full derivative is then promoted to one of the basic variables. Finally, in expression (1901) the underived variable 
ki is used to solve for the Hamiltonian constraint. Note that k^ could also be used to solve for that constraint, but 
this would not obey our third criterion about having a well defined vacuum limit. In this case the perturbations of 
the scalar field would be pure gauge and thus the physical matter degrees of freedom would be encoded in a geometric 
pair. In addition we would not get a master variable that could be consider the generalization of the Zerilli variable. 

In any case, there are other routes that one could follow to construct a polar master gauge invariant. For instance, 
one could also choose for instance a particular background gauge with a fixed (nonzero) value of Hi (or n 3 ) and 
divide by this moment when solving the constraint for (or for tt^). Another alternative could be to relax the first 
condition about the algebraic nature of the transformations. 
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V. EVOLUTION EQUATIONS 

The variable Z (1931) obeys a complicated equation of motion. This section summarizes its differential structure, 
and shows that it is indeed a generalization of the Zerilli equation. In order to simplify the calculations, (t, p = r) 
are chosen as background coordinates. In addition we will take ,5 = 0 which, because of the background evolution 
equation m, implies Hi = 0. By reversing all canonical transformations, it is straightforward to write the variable 
Z in terms of the initial harmonic coefficients (I21H26I) . In particular, in this background gauge the master variable 
takes the following form: 

Z = -—H2 + -^{4:K,r + 2<^,rip - Ulr^G^r) + -^[2K + l{l + 1)G]{12 - 0^(4 + 41(1 + 1) + njr^] - 

CLT’ ^Oj QCLT' 

+ —^ [ 2/(1 + 1 ) + ( 103 ) 

ar^ 

In order to get this expression, it is enough to reverse all transformations presented in the previous section. The 
moments {pi,P 2 tP 3 tP) can be written in terms of time derivatives of their corresponding position variable by making 
use of the linearised evolution equations. Even so, note that all time derivatives disappear from the expression of the 
master variable in the chosen background gauge. 

Let us now present the physical evolution equations. The equations of motion for the master gauge-invariant 
variables are obtained by direct variation of the action (l28)) , 

Pz ^ + ..., (104) 

Zt = -h..., (105) 

P^^^ = M!^f pP + + ..., (106) 

= Mif pP + Z^''^ + Mif + ..., (107) 

where the Roman numerals in the different superscripts represent radial derivatives and the dots stand for terms 
with lower-order radial derivatives. In addition, the subindices of the background-dependent M-coefficients denote 
positions in a matrix. More precisely, would correspond to the slot {i,j) of the matrix that multiplies the fcth 

order radial derivative of the column vector {Pz, Z, P^, Before further analyzing this system of equations, let us 
particularize it to the vacuum case and recover the results obtained in such a scenario by Moncrief [l|. 


A. The vacuum case: the Zerilli equation 


For a vacuum spacetime, all background and perturbative fluid variables {Ha, $, Q^, P 5 } disappear from our problem 
and equations (I106m07'l) are empty. In this case, the background gauge conditions we have chosen at the beginning 
of this section, which implied Hi = 0 , also impose a vanishing value for the background momentum 112 due to the 
constraint (|I3I) . In addition, from background equations (I18I19I) one gets the explicit form of the metric components. 


_ 1 _ 
a 



(108) 


In this way, one can solve equation (jl05l) to write down the gauge-invariant momentum Pz in terms of the time 
derivative of its conjugate variable Z, 


Pz 


l{l + l){2rK + QMf ’*■ 


(109) 


For convenience, we define the rescaled variable, 


Z 

a 2rA -|- 6M ’ 


( 110 ) 


which, inverting the canonical transformations that have been performed in the previous section, can be expressed in 
terms of the initial harmonic coefficients as 

^ = 3M + Ar 


( 111 ) 











11 


Finally, introducing relation (11091) in equation (I104L the Zerilli equation is obtained, 


1 - 


2 M' 


' df^ 


d'^X \ _ 

dr*‘i j 


X = 0, 


( 112 ) 


where we have made use of the tortoise coordinates with r* = r + 2M\w{i^ 

by, 

_l{l + 1) 6 M r2A(A + 2) + 3M(r - M) 

2 ■ ^2 J,3 _j_ 2,MY 


1), and the potential is given 


(113) 


Therefore, the gauge-invariant combination y (11111) reduces to the Zerilli variable when particularized to vacuum. 


B. Simplification of the evolution system 


The highest radial derivative that appears in the evolution equations (llOdHlOTl) is that of the Z variable. More 
precisely, in Eq. (11041) it is of sixth order, whereas in Eqs. (11061) and (11071) it is of fifth order. Nonetheless in Eq. (11051) 
only radial derivatives of the variable Z up to fourth order appear. Thus, this latter equation can be derived with 
respect to r in order to replace those higher-radial derivatives of Z, which appear in other equations, with lower-order 
derivatives of different variables and a time derivative of Z. In fact, it turns out that it is possible to perform a change 
of variables so that the highest radial derivative in the equations is of third order only. Let us define the new variables 


Cl — BiPz -\- B2Z rr + B^Z^r, 
6 = 

Cs = B^Pif, + B^Zr, 

C 4 = Be4> + ByZ^r, 


where the Bi coefficients are given by, 


Bi 

B 2 


B 3 


Bi 


B 5 

Be 


Br 


(rIlaY U2 

a / aD ’ 

a^D ) ’ 


/Usr^V 

2 A2 V a'^D ) 


U 2 




-f r [: 




/(;-f 1)] n2(4>.p)^ 


aV 


n^n; 


2113 


47.2 

Id 


(DUf) U2+r [3a2;(l -f 1) - 8] U 2 UI + -1(1 + 1) [a^ {l^ + I + l) - 3] U 2 


2a fr^UsV U 2 

X \ ) D’ 

/ 2 \ ^ 

f j 112 1 [(1^ -|- ^ -|- 2) -|- 2] $ pn 2 -l- 2 ra^DYie -\- 2 r$_ppn 2 j-, 

—Bi, 


2ar 



Ha. 


The system of equations (I104H1071) is then rewritten in the following way: 

dt 2^ 2^ fc ’ 

j — 1 k—0 

or, by defining the column vector C := (Ci, C 2 , Ca, ^ 4 )^, in matrix notation, 

AT w- ill 

dt ^ ^ dr^ 

fc =0 


(114) 

(115) 

(116) 
(117) 


(118) 

(119) 


( 120 ) 

( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 


(126) 
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The first row of the matrix corresponding to third-order radial derivatives is zero, that is = 0, and thus only up 
to second-order derivatives of appears in the equations. Note that the global differential order of the system of 
equations has been reduced in three (from up to sixth-order radial derivatives to up to just third order) by the above 
transformation (I114mi7l) . 

In order to further simplify this set of equations, one could perform another change of variables that takes the matrix 
coefficients of the higher-order derivatives, in this case A 3 , to its Jordan form. In the axial case such a transformation 
converted a set of two equations of second-order into a new set of an equation of second-order and another with no 
radial derivatives [ 2 ^. 

In the present case, all four eigenvalues of the matrix A 3 are zero and its corresponding Jordan form is, 

/O 1 0 0 
I 0 0 0 0 
~ 0 0 0 1 
Vo 0 0 0 

As it is usual, a similarity transformation, implemented by a matrix S, relates the matrix A 3 with its Jordan form J, 

J = S-^A 3 S. (128) 

With this S matrix at hand, one can define the normal coordinates w as 

w := (129) 

And the equations of motion for these new variables will take the following form: 



dCo d^LO ^ d'^uj 


(130) 


It is easy to see that these n matrices can be written in terms of Ak matrices in combination with derivatives of the 
similarity matrix S in the following way. 


H 2 

= .5-1(3A3^.. + A2^), 

(131) 

Hi 

= S-\3A3S,rr + 2A2S,r+AiS), 

(132) 

Ho 

= S ^ [A3S^rrr T A 2 S^rr “t” AiS J- Ao»S* — ,t') • 

(133) 


Unfortunately, these matrices are quite involve. In particular, and as opposed to what happened in the axial case, 
all the components of the matrix ^2 are non-vanishing and thus second radial derivatives of all variables appear in 
all equations of the system (I130p . Therefore the only advantage of the normal variables w with respect to the initial 
variables ^ lies in the simplicity of the third-order radial derivative terms, which only appear in two of the variables. 
Nonetheless, this advantage might not be so relevant if one considers that our final interest is to obtain the perturbed 
metric in terms of the variables for which we solve. Reconstructing the perturbed metric from the uj variables implies 
another change of variables more than reconstructing it from ^ variables. In addition, matrices are more involve 
than their Ak counterparts for all fc < 2. The lengthy expressions of Ak and prevents us from providing them here 
explicitly. Even so, these matrices are available from the author by request. 


VI. CONCLUSIONS 

In this paper a generalization of the Zerilli master variable for a specific spherical but dynamical background 
spacetime has been presented. In order to factorize and remove the angular dependence from the equations of motion, 
a decomposition on tensor spherical harmonics of the polar part of different perturbative variables has been performed. 
At linearised level each harmonic coefficient, characterized by {l,rn) angular numbers, decouple from the rest due to 
the symmetry of the background. The perturbative problem has then been formulated on a Hamiltonian framework, 
which shows very clearly the dynamical role of each object. As it is well known, the second variation of the Einstein- 
Hilbert action provides an action functional for the perturbative variables. In this way, an effective Hamiltonian can be 
defined for linearised variables, which is given as a linear combination of a physical Hamiltonian and four constraints, 
which are the linearised constraints of general relativity. These constraints are the generators of gauge transformations. 
Thus, in order to construct gauge-invariant master variables, one can perform a canonical transformation so that four 
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of the new variables are equal to the constraints. In this way, for the considered background spacetime, one obtains 
two pairs of variables that encode the complete physical information of the problem: one pair corresponds to the polar 
mode of the gravitational wave and the other one to the scalar matter degree of freedom. 

In order to be admissible, three conditions have been requested to this canonical transformation, which has been 
performed as five subsequent transformations to provide a clear view of each step. First, they should be algebraic 
transformations so that they do not involve any integration and can be performed explicitly. Second, they should 
not require dividing by any background object that could vanish. And third, they should have a well defined vacuum 
limit so that we obtain a generalization of the Zerilli variable. This latter in particular implies that the perturbative 
matter degrees of freedom can not be used to solve for the perturbative constraints, since such a transformation would 
not be well defined in the limit that the perturbations of the scalar field vanish. The full transformation that has 
been proposed here completely fulfills the first and third criteria, but it is unclear whether it also satisfies the second 
one. In the last (fifth) transformation it turned out to be necessary to divide by the background coefficient D defined 
in (EH). This background coefficient is positive definite in the Schwarzschild case (11021) . but it is difficult to assert 
something about it in the general dynamical case. 

In principle the master variable that has been found is not unique, since there is apparently much freedom in 
the canonical transformations that one could perform. Nonetheless, the three imposed criteria reduce considerably 
this freedom. In particular note that, once the canonical transformation (I43H52I) is performed, these criteria almost 
completely single out the subsequent transformations. More specifically, the momenta tts and 7r4 are the only variables 
that can be used to solve for the constraints (IMl) and (15^ respectively. This leads to the form (1551) of the linearised 
Hamiltonian constraint. In that expression none of the variables can be used to solve the constraint algebraically. 
Thus, next parametrized transformation (I68H77I) is performed in order to concentrate all derivatives in a unique full 
derivative. The term inside this full derivative is then promoted to one of the basic variables. Finally, in expression 
EOD the variable fci, which appears with no derivatives, is used to solve for the Hamiltonian constraint. Note that 
could also be used to solve for that constraint, but then this transformation would not have a well defined limit when 
the perturbative matter degrees of freedom vanish and would then violate the third criterion. 

Finally, the evolution equations obeyed by the master variables have been obtained. The differential order, in radial 
derivatives, was initially seven. Nonetheless, it is possible to redefine new master variables so that the highest radial 
derivative is of third order. Hence contrary to the Zerilli variable, which fulfills a wave equation, these master variables 
obey equations of higher (radial) order. This could be surprising since the equation of motion of a linear combination of 
objects, which obey second-order differential equations, is obviously of second-order in terms of the objects themselves. 
Nonetheless, the fact that these objects are expressed in terms of the combination itself (and its radial derivatives) 
can increase the differential order of the equation. This is exactly what happens in this case. Unfortunately the 
obtained equations are quite involved and it is not clear if they could be of practical use. Nevertheless, the fact that 
the mentioned canonical transformations can be performed algebraically and completely decouple the gauge from the 
physical degrees of freedom, turns out to be a relevant result in itself, which could pave the way for the construction 
of a polar master variable for any dynamical spherically symmetric background. 
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